The Roberge- Weiss endpoint in Nf = 2 QCD. 



O 

00 
-f— > 

cd 

Oh' 

> 
in 

in 



X 



Claudio Bonati 1 , Guido Cossu 2 , Massimo D'Elia 3 and Francesco Sanfilippo 4 
1 Dipartimento di Fisica, Universita di Pisa and INFN, 
Sezione di Pisa, Largo Pontecorvo 3, 56127 Pisa, Italy 
2 Theory Center, IPNS, High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801, Japan 
3 Dipartimento di Fisica, Universita di Genova and INFN, 
Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy 
4 Dipartimento di Fisica, Universita di Roma "La Sapienza" and INFN, 
Sezione di Roma, Piazzale A. Moro 5, 00185 Roma, Italy 
(Dated: July 11, 2011) 

We present the results of extensive simulations regarding the critical behavior at the endpoint of 
the Roberge- Weiss transition for Nf = 2 QCD. We confirm early evidence, presented in Ref. [J], 
according to which the Roberge- Weiss endpoint is first order in the limit of large or small quark 
masses, and second order for intermediate masses. A systematic study of the transition strength as 
a function of the quark mass in the first order regions, permits us to estimate the tricritical values 
of the quark mass separating the second order region from the first order ones. 

PACS numbers: 11. 15. Ha, 64.60.Bd, 12.38.Aw 



I. INTRODUCTION 

A full understanding of the QCD phase diagram at fi- 
nite temperature T and baryon chemical potential /is is 
one of the main unreached goals within the Standard 
Model of Particle Physics. Various questions remain 
open, which are of fundamental importance both the- 
oretically and phenomenologically, for astrophysics and 
heavy ion collisions, like the existence and location of a 
possible critical endpoint in the T — (1b plane, accessible 
to experiments. 

Lattice QCD simulations, which are in principle the 
ideal tool for a full non-perturbative investigation of the 
phase diagram, are unfortunately hindered at [ib ^ by 
the complex nature of the path integral measure (sign 
problem). Among other approximate methods, a way 
to partially overcome the sign problem is to consider a 
purely imaginary quark chemical potential, \x q = /is/3 = 
i/ij: numerical simulations are feasible and information 
about real /is can be recovered by analytic continuation 
techniques [34l6j ■ 

Recent literature has pointed out that the phase struc- 
ture at finite T and imaginary chemical potential may be 
important by its own, and teach us something about the 
non-perturbative properties of QCD also at zero or small 
real /xs [Tl [TT1. [l7H20l | . Such phase structure is character- 
ized by a periodicity of the partition function 



Z(T, Mj )=Tr 



(1) 



in the angular variable 9 = fij /T, which can be viewed, in 
the path integral representation of the partition function, 
as a phase rotation of fermion boundary conditions in 
the Euclidean temporal direction. It can be shown [21] 
that the period in 6 is 2ir/N c , where N c is the number 
of colors. Such periodicity is smoothly realized in the 
low temperature, confined phase, as expected from the 
fact that only uncolored states, with N q multiple of N c , 
contribute to the system dynamics. 



The situation is different in the high temperature 
phase, as expected from the fact that also colored states 
appear. Indeed, as can be explicitly verified by pertur- 
bative computations 
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the periodicity is realized in 
a non-analytic way: the system goes through first order 
lines, known as Roberge- Weiss (RW) transitions, when 9 
crosses some fixed values, 9k = (2k + l)ir/N c , where k 
is an integer. For such values of 9 the system possesses 
an exact Zi symmetry, which is spontaneously broken 
for T > Trw and unbroken for T < Trw : therefore at 
T = Trw, which is in fact the endpoint of the RW lines, 
a genuine finite T phase transition takes place for all val- 
ues of the quark masses. Such transition coincides with 
the phase transition at which charge symmetry is sponta- 
neously broken when a spatial dimension is compactificd 
below a given critical size (see e.g. Refs. [22T - [26j for early 
lattice studies of such transition, which has been inves- 
tigated in the context of orientifold planar equivalence 

EE Hi). 

The endpoint of the RW lines has been considered by 
recent literature [J dH Q3, HI H, for its possible 
influence on the critical properties and on the phase dia- 
gram of QCD. The endpoint can be second order in the 
3D Ising universality class, or first order; in the latter 
case it is actually a triple point, from which two further 
first order lines depart. 

In Ref. [l| first evidence has been presented showing 
that, for QCD with two degenerate flavors (Nf = 2), 
the endpoint is first order in the limit of small quark 
masses and second order for intermediate masses; first 
order comes back in the high quark mass regime, where 
the system reaches its quenched limit. In the same paper 
it has been pointed out that, when the endpoint is first 
order (triple point), one of the further first order lines 
departing from it can be identified with (part of) the 
continuation of the critical line to imaginary chemical 
potential, thus explaining early evidence [3, [1] that the 
latter meets the RW line right on its endpoint. A further 
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conjecture, put forward in Ref. lj, has been that the 
nature of the transition at \i = as a function of the 
quark mass spectrum (which is summarized in the so- 
called Columbia plot) is regulated by the physics of the 
RW endpoint itself, i.e. that the fi = transition is first 
order only when the first order line departing from the 
RW triple point reaches the \i = axis. 

Recently the numerical study of the RW endpoint has 
been extended to Nf — 3 QCD [19(, confirming also for 
this case the presence of a first order transition for small 
and high quark masses, with a second order region for 
intermediate masses. Moreover, the authors of Ref. (l9j 
have suggested that the tricritical behaviour which is 
present at the two tricritical masses, separating the sec- 
ond order from the first order regions, may shape the 
critical line also for real values of the chemical potential, 
implying a weakening of the transition with real chemical 
potentials which was suggested also by earlier works [3(| ■ 

All the results and conjectures above claim for a more 
systematic study of the phase diagram in the T — fj,j 
plane, which is perfectly feasible with present simulation 
algorithms. The aim of the present work is to move a step 
in this direction, by extending in a substantial way the 
original results presented in Ref. [l| for Nf = 2 QCD. 
In particular we will present results about the critical 
behavior at the RW endpoint for a large set of quark 
masses, confirming the results of Ref. Ill] and giving an 
estimate for the two tricritical masses, mti and m t 2 > 
TOti, separating the first order regions from the second 
order one. 

Our first instrument to discern the critical behavior 
around the RW endpoint is the finite size scaling of vari- 
ous susceptibilities. However, an accurate determination 
of the critical properties around the tricritical point may 
be a non-trivial task. Much can be learned in this direc- 
tion by the study of simpler statistical systems, like the 
3D 3-state Potts model in presence of a negative mag- 
netic field h [3l], H3] , which shares some of the properties 
of QCD along the RW lines, i.e. the presence of a resid- 
ual Z2 symmetry which gets spontaneously broken at a 
critical temperature. In that model the transition is first 
order for small values of \h\ and second order for large 
values of \h\, with a tricritical value of the field, /itric, sep- 
arating the two regimes 1 . As shown in Ref. [3l| . discern- 
ing the correct universality class close to h tl -i c is difficult 
since, at a given distance from /itrio tricritical scaling 
will mask the correct critical indexes up to a given lat- 
tice size I max , which is regulated by tricritical crossover 
exponents. A similar phenomenon is expected around 
TOti and m t 2- Following Ref. 31], an alternative strategy 
will be to determine parameters which fix the strength of 
the first order transition for m < m t \ or m > ra t 2, like 



1 In the Potts model, of course, one does not observe the re- 
strengthening of the transition (hence a second tricritical point), 
which is present for QCD at low masses and which is likely caused 
by the interplay with chiral degrees of freedom. 



the latent heat or the gap of the order parameter, and 
extrapolate the values of m at which such parameters 
vanish, i.e. the first order transition disappears. 

Our results have been obtained using standard rooted 
staggered fermions on lattices with N t = 4. The paper 
is organized as follows: in Sec. II we give more details 
about the discretized version of QCD under investigation 
and about the observables and the strategy used for the 
study of the critical behaviour; in Sec. Ill we present our 
numerical results and finally, in Sec IV, we discuss our 
conclusions and perspectives. 



II. NUMERICAL SETUP 



We shall consider the partition function of Nf = 2 
QCD in presence of an imaginary chemical potential and 
in the standard staggered discretization of dynamical 
fermions, 



Z(T,6) = JvUe~ SG[u] (det M[U,6]) 



1/2 



(2) 



where 9 = fii/T, Sg is the pure gauge plaquette action 
and AI is the fermion matrix 



Mi. 



1 3 

amS id + - ^2 m,u \ Ui,vSi,j-t> - U}_^ v Si,j+pJ 

+ //..;(«"""/•,.,*.., ; < ■'""I: ,,')„.,) ■ ( 3 ) 

Here i and j refer to lattice sites, v is a unit vector on 
the lattice, r\i. v are the staggered phases, a is the lattice 
spacing and m is the bare quark mass. 

RW transitions take place for 6 = (2k + l)n/3. We 
shall consider in particular the case 9 = ir: for this value 
the residual Z2 symmetry, which is spontaneously broken 
at Trwj corresponds to charge conjugation, hence the 
imaginary part of the Polyakov loop or, alternatively, the 
imaginary part of the baryon number can be taken as 
possible order parameters; as in Ref. [l| , we shall consider 
the former. In the following L will stand for the spatially 
averaged Polyakov loop trace (normalized by N c ), hence 
Im(L) is the order parameter. 

The order parameter susceptibility is defined by 



X = Li (<Im(L)^)-<|Im(L)|>^ 



(4) 



where L s is the spatial size in lattice units, and is ex- 
pected to scale, around the transition, as follows: 



X = li' v mV v ) 



(5) 



where t = (T — Trw)/^rw is the reduced temperature. 
That means that the quantities x/Ll^", measured on dif- 
ferent lattice sizes, should fall on the same curve when 
plotted against tL 1 J v ' . 

Another relevant quantity is the specific heat C of the 
system, which is instead expected to scale as 



C = C + L a s '» fc(tLl /u ), 



(6) 
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3D Ising 


0.6301(4) 


1.2372(5) 


0.110(1) 


~ 1.963 


~ 0.175 


Tricritical 


1/2 


1 


1/2 


2 


1 


1 st Order 


1/3 


1 


1 


3 


3 



TABLE I: Critical exponents relevant to our analysis. 

where Co is a regular contribution. The values of the crit- 
ical indexes a, 7 and v which are relevant to our analysis 
are listed in Table|T](see e.g. Refs. H^HH), together with 
the values they take for the different critical behaviors 
which may take place in our system, i.e. first order, sec- 
ond order in the universality class of the 3D Ising model, 
and tricritical mean field. 

A careful verification of Eqs. ((5$ and ((BJ, as well as 
of similar relations giving the finite size scaling behav- 
ior of other relevant quantities, gives information about 
critical indexes, hence about the universality class of the 
transition. A more direct way, in the case of a first order 
transition, is to verify the existence, in the thermody- 
namical limit, of finite gaps in the order parameter or in 
the internal energy (latent heat), which may be visible 
by looking at double peak distributions of physical ob- 
servables around the transition, or by studying the large 
volume limit of some cumulants. 
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FIG. 2: Reweighted distribution of the real part of the 
Polyakov loop at the pseudo-critical point for am = 1.5 and 
various lattice sizes. 
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FIG. 1: Monte-Carlo histories of the real and imaginary part 
of the Polyakov loop for a /3 value (5.328) around the critical 
point and am = 0.0175 on a 16 3 x 4 lattice. 



An example is the Binder-Challa-Landau cumu- 
lant [35| of the energy, which is defined as B4 = 1 — 
(E 4 )/(3(E 2 ) 2 ). It can be shown (see e.g. [36]) that near 
a transition develops minima whose depth scales as 



BA 
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+ <W) 



(7) 
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FIG. 3: Reweighted distribution of the real part of the 
Polyakov loop at the pseudo-critical point for am = 1. and 
various lattice sizes. 



of B 



-EL). In particular, the thermodynamical limit 
, is less than 2/3 if and only if a latent heat 
is present. To simplify our analysis we have considered 
the average plaquette (sum of the spatial and temporal 
plaquettes) in place of the internal energy, since it is a 
quantity which can be measured much more easily and, 
like the internal energy, is even under the Zi symmetry 
which gets broken at the RW endpoint. To simplify the 
notation, in the following we will use the shorthand 

2 
3 



B 



B 



4 \ min 



(8) 



and from Eq. ([7]) it follows that B oc A|;, where in our 
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FIG. 4: Reweighted distribution of the real part of the 
Polyakov loop at the pseudo-critical point for am = 0.5 and 
various lattice sizes. 
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FIG. 5: Reweighted distribution of the plaquette (average of 
spatial and temporal) at the pseudo-critical point for am — 
0.005 and various lattice sizes. 



case by A^ we actually mean the gap at the transition 
in the average plaquette. 

A different, but analogous quantity is the gap of the 
order parameter, A, which can be extracted by looking 
at the scaling of the maximum of its susceptibility, 
and using the relation, valid in the large volume limit for 
a first order transition, 



Xn 



const. 



Li 



A 



(9) 



Both A e and A are expected to vanish as we approach 
a tricritical mass mtric from the first order side. In partic- 
ular, the leading order expected behavior is the following 
(see 37 1 or [38[ for a brief summary) 



and 



A E CX y/h- /itric 



A OC y/\(h - /l tr ic) l0g{h - htri c )\ 



(10) 



(11) 



where we have indicated generically by h the relevant pa- 
rameter driving the change from first to second order. It 
is clear that h is a function of the quark mass and that 
close enough to the tricritical point one can always set 
h — ft-tric ~ m — Titric! however, appropriate choices of h 
can improve the region around the tricritical mass where 
Eqs. pU|) and (fTTj) hold. Our choice will be h ~ m in the 
low mass region and h ~ 1/m in the high mass region. It 
is interesting to notice that Eq. (jlip may seem ambigu- 
ous, since a multiplicative redefinition h — > const, x h 
changes the functional dependence; however, as long as 
(h — /ltric) <C 1, the change is subleading and Eq. (TIT]) 
still gives the dominant contribution. 

Close to the tricritical points it can be particularly 
difficult to discern the correct critical behavior taking 



place in the thermodynamical limit. Indeed, while first 
order/3D Ising scaling are expected to take place for a 
continuous range of values of m and exact tricritical scal- 
ing only for specific values m = m tr ; c , what really hap- 
pens is that tricritical scaling regulates a neighborhood 
of m tr ic, whose size goes to zero as L s — > oo according 
to critical indexes known as crossover exponents (see e.g. 
[34l [39I I40I]). Indeed, the true critical behavior of the 
system can be seen only for \t\ < \h — /itricl 1 ^, where t is 
the reduced temperature and <p is the crossover exponent, 
which is by definition <j> — y^ /y t (yt and yt are the renor- 
malization group eigenvalues of the relevant variables t 
and h — h tI i C ), in particular (j) = 1/2 in our case [37| . 
Putting the question the other way around, on a finite 
lattice of typical size L s , \t\ can be traded for L s and 
the previous condition becomes L s > \h — /itricl"^; in 
particular, according to the known tricritical indexes in 
Table HI one expects tricritical behavior to dominate and 
mask the correct thermodynamical limit up to a critical 
size 



L r 



A\h- h t 



(12) 



where A is some unknown constant. Such a behavior has 
been studied and verified quantitatively in Ref. (3ll | in the 
case of the 3D 3-state Potts model in a negative external 
field, which shares part of the symmetries studied in the 
present work. 

The difficulties in discerning the correct critical behav- 
ior around m tr ic may result in a difficult determination 
of the tricritical mass itself. For this reason we have fol- 
lowed the strategy adopted in Ref. [3l|, i.e. to determine 
the cumulant of the plaquette B and the gap of the order 
parameter A 2 for values of m where a first order tran- 
sition is present, and then to determine m tr i c by fitting 
data with the expected behaviors reported in Eqs. (TTU)) 
and (fTT|) . 
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FIG. 6: As in Fig. for am = 0.01. 
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FIG. 8: Monte-Carlo histories of the Polyakov loop (abso- 
lute value) and of the chiral condensate for a /3 value (5.314) 
around the critical point and am = 0.01, on a 16 3 x 4 lattice. 
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FIG. 7: As in Fig. [5] for am = 0.075. 



am 


B 


A 2 /4 


0.005 


2.15(10) x 10" 4 


9.60(20) x 10" 3 


0.010 


1.54(7) x 10~ 4 


8.04(26) x 10~ 3 


0.0175 


1.01(8) x 10" 4 


6.40(40) x 10" 3 


0.025 


0.69(4) x 10" 4 


5.54(24) x 10 -3 


0.030 


0.48(7) x 10" 4 


4.60(50) x 10" 3 


0.035 


0.32(6) x 10~ 4 


3.60(40) x 10" 3 


1.00 


0.38(4) x 10" 5 


2.59(13) x 10" 3 


1.25 


0.58(7) x 10~ 5 


4.16(36) x 10" 3 


1.50 


0.66(7) x 10" 5 


4.32(24) x 10" 3 


2.00 


0.89(7) x 10~ 5 


5.20(20) x 10" 3 



TABLE II: Estimated values for the thermodynamical limit 
of B and A 2 /4 for values of the quark mass where a first 
transition takes place. 



With the aim of determining the tricritical masses 
TOti and 777.42 present in the low and high mass regions 
respectively, we have studied the critical behav- 
ior of the system for various quark masses, am = 
0.005, 0.01, 0.0175, 0.025, 0.03, 0.075, 0.2, 0.5, 1., 1.25, 1.5 
and 2.0. For each quark mass we have made simulations 
on lattices with N t = 4 and different spatial sizes L s , 
reaching up to L s = 40 when necessary to correctly 
discriminate the critical behavior. Numerical simulations 
have been performed using the standard Rational Hybrid 
Monte-Carlo algorithm 41]. Collected statistics have 
been typically of the order of 10 5 trajectories around the 
critical f3 and for each value of L s . 

Apart from results obtained for am = 0.025 and 
am = 0.075, which were already partially reported in 
Ref . [l| , most numerical simulations have been performed 
on two GPU farms located in Pisa and Genoa and pro- 
vided by INFN, consisting of a total of 8 S1070 (32 
C1060) NVIDIA GPUs. The numerical code, which runs 



almost entirely on the GPUs, has been described in detail 
in Ref. H|. 

III. NUMERICAL RESULTS 

The presence of a first order RW endpoint, i.e. of a 
triple point at the end of the RW lines, has clear sig- 
natures in the Monte-Carlo (MC) histories and in the 
probability distributions of the order parameter and of 
other quantities. In Fig. Q] we show the MC histories 
of the real and imaginary part of the Polyakov loop for 
am = 0.0175, where the endpoint is first order, and a j3 
value around the transition. Metastabilities are clearly 
detectable, with Im(L), the order parameter, taking three 
distinct possible values, one in the unbroken and two in 
the broken Z2 phase. Re(L), which is Zi even, takes in- 
stead only two distinct values corresponding to the bro- 
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FIG. 9: Scaling of the reweighted susceptibility of the imaginary part of the Polyakov loop according to first order critical 
indexes for am — 0.0175 (up-left), am — 0.03 (up-right), am — 1.5 (down-left) and am = 1.0 (down-right). 




FIG. 10: Scaling of \ f° r am = 0.075 according to 3D Ising critical indexes (left) and to tricritical mean field indexes (right). 



ken and unbroken phase. 

In Figs. [2l [3] and [4] we show the reweighted distribu- 
tion of Re(i), at the pseudocritical values of /3 taking 
place on the different lattice sizes, for three values of 
am in the heavy quark region, am = 1.5, 1.0 and 0.5 re- 



spectively. For am — 1.5 and am = 1.0 a double peak 
distribution clearly develops and deepens as L s — > oo, 
indicating a first order transition, even if in the latter 
case one has to reach L s = 40 to clarify the behavior, 
indicating that in this case the first order transition is 
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FIG. 11: Maximum of the susceptibility of the real part of 
the Polyakov loop as a function of the lattice size L s and for 
am — 0.2, together with a linear fit including sizes L a < 32. 



13: Binder-Challa-Landau cumulant of the plaquette 
definition in Eq. j5}) as a function of the lattice size 
0.025 and am = 0.075. In the first case a func- 

1 well all data with 

For am — 0.075, in- 



FIG 

(see 
for am 

tion B = a + b/Li + c/L° a describes 
a = 0.69(4) x 10 -4 and xV d -°- f - = 0.13 



stead, data with L s > 8 are well described (x 2 /d.o.f. = 0.69) 
by a dependence B — aL h s (b — 0.62(2)) which gives B — in 
the thermodynamical limit. 
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FIG. 12: Maximum of the susceptibility of the order param- 
eter, x, as a function of the lattice size L s for am = 0.025, 
together with a cubic fit \ = const. + A 2 including all 
sizes (x 2 / d -°-f- = 0.89). 



weaker. For am = 0.5, instead, the distribution stays 
single peaked for all explored volumes, suggesting that 
the endpoint may be second order in this case: this hy- 
pothesis is indeed consistent with the determination of 
a.iTit2 presented later. 

Similar considerations can be made for the light mass 
region. In Figs. [SJ H] and [7] we show the reweighted pla- 
quette distributions at the pseudocritical couplings for 
am = 0.005,0.01 and 0.075 respectively. Double peak 
distributions are present for the two lower masses, with 
the first order being clearly stronger for am — 0.005. 



For am = 0.075 instead, as already shown in Ref. [l|, 
the distribution stays single peaked, suggesting that the 
endpoint is second order in this case: this is consistent 
with our determination of am t \ (see later). 

It is interesting to notice that, when the transition is 
first order, a gap develops also in other quantities, includ- 
ing the chiral condensate, as visible from Fig. [HJ where we 
show the MC histories of the chiral condensate and of the 
Polyakov loop around the RW endpoint. That suggests 
that, as for the usual thermal transition at fj, = 0, a strict 
correlation between deconfinement and chiral symmetry 
restoration may be present also at the RW endpoint. 

These results already fully confirm the outcome of 
Ref. [l| : the RW endpoint is first order in the chiral limit 
and weakens as the quark mass is increased, till an in- 
termediate mass region is reached where the transition is 
second order; it is first order again in the high quark mass 
limit, where it weakens as the quark mass is decreased. 
Last result is in some sense trivial since, as already dis- 
cussed in Ref. [y, it is expected from the fact that the 
SU (3) pure gauge transition is first order. 

Further confirmations come from looking at the finite 
size scaling of the susceptibility of the order parameter, 
X, which is shown in Fig.[9]for am = 0.0175, 0.03, 1.5 and 
1. The first order scaling ansatz, Eq. (jS"|) . is always veri- 
fied for the largest volumes available. However, typically 
one has to go beyond some critical size before seeing the 
correct asymptotic critical behavior, and this critical size 
increases as the transition weakens, i.e. as we approach 
the tricritical points. For instance, at am = 1 first order 
scaling sets in only for L s > 32. 
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FIG. 14: Binder-Challa-Landau cumulant of the plaque- 
tte, extrapolated to the thermodynamical limit, and A 2 /4 
for small quark masses where a first order transition is 
present. We include the result from a linear fit B^ — 
b (arriti — am), giving the value of the tricritical mass am t i = 
0.0428(24) and xVd.o.f. = 0.13 (we have included quark 
masses am > 0.0175), and from a fit to Eq. (fTTI) . A 2 /4 = 
c (am t i — am) \og(am t i — am), giving am t \ = 0.0477(23) and 
X 2 /d.o.f. = 0.37 (we have included quark masses am > 0.01). 
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FIG. 15: Binder-Challa-Landau cumulant of the plaquette, 
extrapolated to the thermodynamical limit, and A 2 /4 for high 
quark masses where a first order transition is present. We in- 
clude the result from linear fits = b (l/(amt2) — l/(am)), 
giving am t2 = 0.71(4) (x7 d -°- f - = 1-09), and A 2 /4 = 
c (l/(am t2 ) — l/(am)), giving am t2 = 0.67(3) (x 2 /d.o.f. = 
1.0). All masses have been included in the fit in both cases. 



Similar considerations apply to the second order re- 
gion. On the left-hand side of Fig. [TU1 which is taken 
from Ref. |l|, we show the finite size scaling of x f° r 
am = 0.075 according to 3D Ising critical indexes: scal- 
ing is fair for the heights of the peaks and less fair for the 
widths. On the contrary, we realize that tricritical mean 
field indexes perform much better, as apparent from the 



right-hand side of Fig. [10] (notice from Table Q] that j/v, 
regulating the height of the peaks, is practically the same 
for 3D Ising and tricritical mean field, while 1/V, which 
regulates the widths of the peaks, is different). That does 
not mean, of course, that am — 0.075 is exactly equal to 
one of the two tricritical masses, but rather that it is close 
enough to one of them so that a fake tricritical scaling 
masks the correct asymptotic scaling at least for sizes up 
to L s — 32. However, we do not know neither how close 
we are to the tricritical mass, nor how large we have to 
go with L s to reach the thermodynamical limit, since we 
have no apriori knowledge of the prefactor appearing in 
Eq. (12]). 

A quantity which is well suited for discerning 3D Ising 
from tricritical behavior is the specific heat C. Indeed the 
coefficient a/v, which regulates the scaling of the height 
of the singular part of C (see Eq. (J6])), changes apprecia- 
bly when going from tricritical to 3D Ising critical behav- 
ior (see Table fl} , hence deviations from tricritical scaling 
are expected to appear first in such quantity. A direct 
measure of the specific heat of the system is not an easy 
task, however the susceptibility of any quantity, sharing 
the same transformation properties of the energy under 
the relevant Z2 symmetry, is expected to scale in the same 
way: examples are given by the plaquette or by the real 
part of the Polyakov loop, which are both Z2 even. In 
Fig. [11] we show the susceptibility of the real part of the 
Polyakov loop as a function of L s for am = 0.2, which we 
expect to be in the 3D Ising region. It is apparent that 
data follow a linear behavior (i.e. a/v = 1), with devia- 
tions visible only for L s > 32 and going in the direction 
of a smaller value of a/v (as expected for 3D Ising); in 
particular in the figure we have plotted the result from a 
linear fit to data up to L s = 28. 

Therefore, in order to get a more reliable determination 
of the tricritical masses, we follow the strategy described 
in Sec. II and proceed to a determination of the gap of the 
order parameter and of the plaquette as a function of the 
quark mass in the first order regions. In Fig. we plot 
the maxima of the order parameter susceptibility, \i as a 
function of L s , for am = 0.025, together with a fit to the 
asymptotic expected behavior, Eq. ©, from which we 
extract A 2 /4, The same procedure has been repeated for 
all quark masses where a first order transition is present. 
In Fig. [T31 instead, we plot the Binder-Challa-Landau 
cumulant of the plaquette, B (see Eq. ((5|)), as a function 
of 1/V for am = 0.025 and am = 0.075: in the first case 
the cumulant extrapolates to a non-zero value as V — > 00, 
with both linear and quadratic corrections in 1/V clearly 
visible, while in the second case data are well described 
by a power law and B — as V —> 00, indicating the 
absence of a gap in the plaquette. 

In Table |TT] we summarize all determinations obtained 
for B and A 2 /4. From such values we can try to deter- 
mine the tricritical masses as the points where B and A 
vanish, fitting data to the expected behaviors shown in 
Eqs. (|TUJ) and ([IT]). In Fig. H31 we show the results of 
such fits in the low mass region for B and A 2 /4, respec- 
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tively. We obtain am a = 0.0428(24) from B. Instead, 
from A 2 /4, we get am tl = 0.0477(23) if we fix h = m 
in Eq. (jlip . however in this case one should take into 
account also the systematic uncertainty related to a pos- 
sible multiplicative redefinition, h — Ahm. In order to 
further check that our results for B and A 2 /4 can indeed 
be described in terms of a common tricritical mass, we 
have also performed a combined fit to all data obtained 
in the low mass region according to 

B = b (am — amti) 
A 2 /4 = c (am — amti) log(Ah(am — amti)) ; (13) 

including directly, in this case, the possible multiplicative 
redefinition h = Ahm among the fit parameters. The 
best fit gives b = -387(46), c = 0.17(6), A h = -9(5) and 
m t i = 0.043(2), with a x 2 /d.o.f. = 0.3/4: the hypothesis 
is therefore well verified, but we cannot trust the uncer- 
tainties on the parameters deriving by this best fit, since 
data for B and A 2 /4 are correlated; notice also that the 
multiplicative constant Ah is poorly determined. Staying 
conservative with the error estimate, we take as our final 
determination m t i = 0.043(5). 

In Fig. [15] we show instead the same kind of fits for the 
high mass region: in this case we have used l/(am) as 
the relevant variable h, as explained in Sec. II. We obtain 
am t 2 — 0.71(4) from B. Instead, regarding A 2 , we notice 
that (h — /itric) is 0(1) and it makes no sense to look for 
logarithmic corrections (see Eq. (fTTjl ): a simple linear fit 
for A 2 (see Fig. [T5|) gives am t 2 = 0.67(3). However, also 
in this case we can redefine h = Ah/m and try again a 
combined fit according to 



B = 



H-L--L- 

am amti 



(14) 



A 2 /4 



1 



1 



81)1(2 



log [A 



1 

am 



am t 2 



leading to am t2 = 0.72(5) and Ah ~ 10~ 2 , with 
X 2 /d.o.f. = 2.2/4. Also in this case one should take into 
account correlations among data for B and A 2 /4, hence 
we prefer to stay conservative in our error estimate and 
state am t2 = 0.72(8). 

We notice that both determinations, am t i — 0.043(5) 
and amt2 — 0.72(8), are consistent with the fact that 
the quark masses for which no metastabilities and double 
peak distributions are observed (am = 0.075, 0.2, 0.5) are 
within the second order region. 



IV. CONCLUSIONS AND PERSPECTIVES 

We have confirmed the outcome of Ref. [l| regarding 
the order of the endpoint of the RW transition for Nf = 2 
QCD: a first order endpoint (triple point) is present both 
in the low mass and in the high mass limit; the endpoint 
is second order for intermediate quark masses, which are 
separated from the first order regions by two distinct tri- 
critical masses. Following an investigation performed in 



ti 



m r 



'a 



FIG. 16: Sketch of the phase diagram in the T-m q plane 
which summarizes our results: in Nf = 2 QCD the endpoint 
of the Roberge- Weiss transition is first order close to the chiral 
and to the quenched limit and second order for intermediate 
masses. A conservative estimate for the two tricritical masses 
separating the second order region from the first order ones, 
for the lattice discretization adopted in the present work, is 
amti = 0.043(5) and arn t2 = 0.72(8). 



Ref. [31| for the 3D 3-state Potts model in a negative 
external field, which shares part of the same symmetries 
studied in the present work, we have performed a careful 
study of some parameters directly linked to the strength 
of the first order transition, in particular the Binder- 
Challa-Landau cumulant of the plaquette and the gap 
of the order parameter; that has permitted to obtain 
independent and consistent determinations of the two 
tricritical masses. Staying conservative with error esti- 
mates, we state as our final result am t i — 0.043(5) and 
am t 2 — 0.72(8). Such results are summarized in Fig. 1161 
where we sketch a phase diagram in the T-m q plane. 

The value of amti corresponds to a pion mass of the 
order of 400 MeV, hence we conclude that for physical 
quark masses the RW endpoint should be well inside the 
first order region. It is therefore of primary importance 
to explore what is the fate of the further first order lines 
departing from the triple point. One of them, in partic- 
ular, may reach the zero density axis or have a critical 
endpoint arbitrarily close to it, which could have great 
influence on the physics of strongly interacting matter 
right above the deconfinement transition. The question 
is also strictly connected to the problem of the order of 
the chiral transition for Nf = 2 |43|, [44| . 

Another important issue is of course to extend our in- 
vestigation to Nf 7^ 2 and confirm the conjecture that 
the nature of the transition at fi = may be regulated 
by the physics of the RW endpoint [1J , i.e. that the fi = 
transition is first order only when the first order line de- 
parting from the RW triple point reaches the = axis, 
and that tricritical scaling indeed shapes the chiral criti- 
cal surface [HI . 

All these investigations will require extensive numeri- 
cal simulations, which are however perfectly feasible since 
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they involve an imaginary chemical potential. Part of this 
program is progress. 

We stress that our present results are valid for the stan- 
dard rooted staggered discretization of the theory and 
for lattices with N t = 4, corresponding to a lattice spac- 
ing of about 0.3 fm. A key issue is then also to verify 
that the main features of the phase diagram remain un- 
changed when changing discretization and/or approach- 
ing the continuum limit. The two tricritical masses could 



still be present, but the first order regions could in prin- 
ciple extend or shrink in a significant way. 
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